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\^ ' A Hamilton-Jacobi formulation of the Lyapunov spectrum and KS entropy 

is developed. It is numerically efficient and reveals a close relation between 
the KS invariant and the classical action. This formulation is extended to 
!>■ , the quantum domain using the Madelung-Bohm orbits associated with the 

Schroedinger equation. The resulting quantum KS invariant for a given orbit 



00 



O ' equals the mean decay rate of the probability density along the orbit, while 
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its ensemble average measures the mean growth rate of configuration-space 
information for the quantum system. 
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Chaos is a long-term instability exhibited by a dynamical system. It is quantitatively 
measured by the Lyapunov spectrum of characteristic exponents, which represents the prin- 
cipal rates of orbit divergence in phase space, or alternatively by the Kolmogorov- Sinai (KS) 
invariant, which quantifies the rate of information production by the dynamical system [|l|. 
The KS invariant is a fundamental quantity that underlies the information dynamics of a 
broad class of dynamical systems well beyond classical mechanics. The results reported 
here show that it is intimately related to the classical action for Hamiltonian systems and 
to amplitude decay and configuration-space information for quantum systems. While chaos 
is a ubiquitous feature of nonlinear classical dynamics, it is conspicuously absent in finite 
quantum systems as these are known to be dynamically stable 0. Nevertheless, the chaotic 
nature of a given classical Hamiltonian produces certain characteristic features in the dy- 
namical behavior of its quantized version. These features have been studied extensively and 
are commonly referred to as "quantum chaos" 0,§]. They include short-term instabilities 
and diffusive behavior versus dynamical localization, characteristic spectral statistics, wave 
function scarring, and other effects. A quantitatively precise formulation of quantum dy- 
namical instability smoothly straddling the quantum-classical transition, on the other hand, 
is at present lacking and is the subject of this work. 

This Letter presents a general approach to the information dynamics of the quantum- 
classical transition based on the Hamilton- Jacobi (HJ) formalism. The formulation is exact 
and based on phase-space concepts, with the KS invariant playing a central role. The exten- 
sion to the quantum domain is accomplished by means of the orbits introduced by Madelung 
and later rediscovered by Bohm . These quantum orbits are natural extensions of the 
classical phase-space flow to quantum mechanics, and provide the required bridge across 
the transition. The present approach yields several new results, including (i) a powerful 
and numerically efficient method for calculating the KS invariant for classical Hamiltonian 
systems, (ii) the striking finding that the quantum KS invariant for a given Madelung-Bohm 
(MB) orbit is equal to the mean decay rate of the probability density along the orbit, (iii) the 
result that the quantum KS invariant averaged over the ensemble of MB orbits equals the 
mean growth rate of configuration-space information, and (iv) a general and rigorous argu- 
ment for the conjecture that the standard quantum-classical correspondence (or the classical 
limit) necessarily breaks down for classically chaotic Hamiltonians. Several other results, 
including those of numerical simulations for classical and quantum models, are presented 
and discussed. 

Canonical formulation of chaos. Our objective here is a symplectically reduced 
formulation of the Lyapunov spectrum for a Hamiltonian system in terms of the classical 
action; see Ref. [[7|, hereafter referred to as "Paper I," for definitions and notation. Consider 
a classical system of A^ degrees of freedom described by the canonical variables {qi,Pi}, i = 
1, . . . ,N, the Hamiltonian function -ff (q, p, t), and Hamilton's principal function ^(q, t, po), 
with Po denoting the initial momenta [§]. Hamilton's equations can be compactly stated 
in matrix form as ^ = J"V^H{^,t), where $, stands for the 2A^-dimensional phase-space 
vector (q, p). Here JT" is a real, antisymmetric matrix of order 2N with a 2 ® N block 
form {On, In, —In, On), which is simply a hsting of its blocks in the order (11, 12,21,22). 
The tangent dynamics of the system is described by the 2A^ x 2A^, nonsingular matrix 

'^!/(^)^o) = 9^iJ.{'tT^o)/9Cou, the sensitivity matrix, where ^(t,^o) i^ the trajectory that 
starts from ^q at t = 0. 



Our first task is to express T in terms of Hamilton's principal function 5*. Using the 
properties of S, we find the desired result in the block form 

-^ — l*^poq ; ~*^poq *^popo' 

*-'qq*-'poq ! *^poq ~ '^qq'^poq '-'poPO/' \^) 

where (Sp^q) . . = d'^ S / dq.dp^. 0. Note that T2iTu^^ is a symmetric matrix (it equals 5'qq), a 
fact that is a direct consequence of the symplectic symmetry of the underlying Hamiltonian, 
which symmetry implies the condition TJ'T = i7 0. This allows us to transform the 
sensitivity matrix T to an upper-triangular block form F according to F = Q{Q)T, where 
f2(9) = (cos 9, — sin 6, sin 6, cos 6) and tan(O) = — Sqq [|Tn[. Note that 6 is a real, 
symmetric matrix of order A^, while Q is orthogonal and symplectic. We shall refer to B as 
the symplectic phase matrix for reasons that will become clear in the following. 

The upper-triangular form (Fii,Fi2, On, T22) of F, on the other hand, guarantees that 
Til = ^22, and that the upper (lower) half of the Lyapunov spectrum is obtained from the 
singular values of Fn (F22) as explained in Paper I. In particular, the KS entropy is given 
by El 



fc = limt_oolog[det(Fn)]/t. (2) 

The evolution equation for T, Eq. (^ of Paper I, leads to the set of matrix equations 
tiiTii^ = L^ + IC^i, Q = CQ, L^ + IC^i = 0, (3) 

which can be used to calculate the Lyapunov spectrum efficiently. Here £ = 
{L^,L^,-L^,L^), /C(t,^o) = VgV5^[^(t,^o).^]. and /C^ = niCCl, where L^ and L^ 
are symmetric and antisymmetric matrices respectively. In particular, we note that the 
time average of the trace of the first equation in the set yields the KS invariant, i.e., 
k =< tr(riiF]3^'^) >, where < / >= lim(_,ooi~^/o/(^'^^- Using the above set of equations, 
we find 

a + /Cii + /C120- + 0-/C21 + 0-/C220- = 0, (4) 

fc =< tr[ sm2BH coszBJ >, (5j 

where a{t) is the restriction of Sqq = — tan(0) to a specific orbit and cr(0) = 0. Equa- 
tions (^ and (^, together with the equations of motion for the orbit, constitute a closed 
set from which the KS invariant can be calculated for any Hamiltonian system. The main 
computational burden is in solving the N{N + l)/2 differential equations in Eq. @). As with 
the results of Paper I, the above set of equations is free of exponentially growing quantities 
and does not require any orthonormality maintenance. 

To illustrate the power of these equations, let us consider the standard form H = p^/2 + 
V{q,t), where q and p are A^-dimensional vectors. We then find K-u = VgVgV^(q, t), 
/C12 = /C21 = On, /C22 = In, and 

cr + a" + /Cn = 0, (6) 



A; = 1 < tr(/Cii - 1) sin 26 > . (7) 

Moreover, using o = — tan(0) and Eq. (|^), we can rewrite Eq. (^ in the remarkably simple 
form 

k = < trcr >p.^,, (8) 



where "p. v." stipulates a principal-value evaluation [T^. Since tr(a) equals V^S along the 
orbit, Eq. (H) simply states that the KS invariant equals the time-average of the Laplacian 
of the action along the orbit. 

We will consider two examples here. The first is defined by A^ = 3 and V{q) = 
-|[(9i - Q2f + (g2 - qsf + {q-i - qif] + {qf + 2g| + 3g|). Equations (§) and (0) together 
with Hamilton's equations, 13 in all, were integrated for a duration of 2.2 x 10^ time units 
with the initial value of H set at unity |TB[. The result was k = 0.6126, with the fluctuations 



never exceeding 0.0001 as the simulation was extended by one order of magnitude to the 
endpoint. In fact a 1% result was already achieved at about 10^ time units. 

For a second example, we will derive explicit formulas for an iV- dimensional kicked system 
(which is otherwise free) with V{q,t) = f{q)Y^'^=iS(t — nT). Let (qn,Pn) be the coordinates 
just after the n*^ kick. Then (qn+i, p„+i) = (q„, p„) + [p„T, -V/(qn+i)]. From Eq. (|), on 
the other hand, we find d+o"^ = for t ^ uT and Acr„ = — VV/(qn) across the n^^ kick. The 
corresponding solution is easily found to be a{t) = (t — uT + a~^) , nT < t < {n + 1)T. 
Thus the iteration map for a is 

an+i = (a-' + T)-' - VV/(q„+i), (9) 

where cr(0) = 0. We now use Eqs. (P) and @ to find 



which, together with the iteration maps given above, provides an efficient algorithm for 
calculating the KS invariant for kicked systems. As an illustrative example, we chose /(q) = 
-l[(gi - q2f + (g2 - qsf + (gs - qif] + {qf + ql + ql), with T = 1.0 x lO'^o. The value 
k = 1.5 X 10^ was obtained in about 10^ iterations [Q. 

Note, incidentally, that in case an eigenvalue of a~^ is negative and less than T in 
magnitude, a will encounter a simple pole in the n + P* cycle [other types of singularities 
are ruled out by Eq. (^)]. This behavior is generic for chaotic Hamiltonians (and common 
for all Hamiltonians), since otherwise a bounded a would, by virtue of Eq. (§), lead to A; = 0. 
Correspondingly, the symplectic phase angles (eigenvalues of 9) go through (m — i)7r as t 
goes through the m^^ singularity of that eigenvalue [we take 0(0) = 0]. These poles occur 
when Til becomes singular (as a matrix), and are the analogues of the so-called conjugate 
points which occur when T12 becomes singular ^. Further insight is gained by considering 

these angles for quadratic potentials for which /Cn is a constant, positive-definite matrix 

i_ 1 

(and of course k = 0). The solution to Eq. (||) is found to be a = — /Cn tan(/Ciit). Thus the 
symplectic phases are given by 9i = arctan[ti;j tan(ti;jt)], where {ui} are the characteristic 
frequencies of the system. Clearly, the mth singularity associated with 6i occurs for Uit = 
{m — 1)77 as expected. The conjugate points occur midway between the latter at ujit = mn. 



Extension to Quantum Mechanics. The MB formahsin plH] associates a phase space 



flow with a quantum system by setting ]T5 



V'lx, t) = exp[iS(x, t)/h + R(x, t)], (11) 

and treating S as Hamilton's principal function for an ensemble of orbits with conflguration- 
space density exp(2R). These orbits are deflned according to 

q(t, qo, Po) = p{t, qo, Po) = VS[q(t, qo, Po),t], (12) 

where qo=q(0, qo, Po)- Here ip{'x,t) is a solution of the Schroedinger equation with the 
Hamiltonian p^/2 + l^(x, t). It is readily verifled that the expectation value of any observable 
in the state ilj{'x,t) is given by its average over the ensemble of orbits defined above. In 
particular, the ensemble average of Eq. (|T^) will lead to Ehrenfest's equations. 

The correspondence thus established allows us to define the quantum KS invariant for a 
given orbit as 

k 4|f < V^S >,.„„ (13) 

in strict analogy to Eq. (^). It is worth emphasizing here that the averaging process in 
Eq. (|13D is with respect to the time along the MB orbit to which S is restricted. 

Intuitively, one would expect that orbits neighboring a hypothetical chaotic orbit in the 
ensemble diverge from it on the average, thus causing the orbit density along the chaotic 
orbit to decrease with a mean rate related to k. Remarkably, this expectation is fully realized 
by the Schroedinger equation. To see this, consider the equation of motion for R(x, t) as 
inherited from the Schroedinger equation: c}R/9t + VR-VS = — |V^S. An inspection of this 
equation shows that its characteristic curves are the MB orbits, so that it takes the following 
form along those orbits: 

rfR(q(t,qo,Po),t)M = --V'S[q(t,qo,Po),t]. (14) 

Using this in Eq. ([T3|) , we find k= — 2< dR/dt >p^ , or equivalently 

k= - < rfln|V[q(t,qo,Po),t]|'M >p...- (15) 

This exact result states that the quantum KS invariant for a given orbit is the mean decay 
rate of the probability density along the orbit. We mention in passing that if ip in Eq. (|r^) is 

replaced with its classical limit ips.c. = [det(S'p„q)]2 exp[iS{'x,t)/h], the result is the classical 
KS invariant k 0. 

Let us consider the implications of Eq. ([I5|) for classically chaotic systems, for which 
k ^ while k = 0. For a chaotic orbit in classical dynamics, there is no integral of motion 
in phase-space except possibly energy. For the ensemble of MB orbits, on the other hand, 
the defining equations (^2|) provide A^ such integrals (upon eliminating q) if VS is a long- 
term predictable, or computable, function. But VS is computable, since the Schroedinger 
equation for finite-dimensional, bounded systems is known to be free of sensitivity to initial 
conditions 0. The inevitable conclusion is that while VS (classical) is not computable for 
chaotic systems, VS (quantum mechanical) is, regardless of how small h may be. Clearly, 



the statement that VS reduces to VS as ^ — ^ 0, often referred to as the classical limit, 
cannot hold for chaotic Hamiltonians. Recalhng that regular Hamiltonians are very special 
while chaotic ones are by far the more typical class, we recognize the very limited validity 
of the classical limit. 

The above considerations confirm the view that the absence of sensitivity to initial con- 
ditions in the quantum dynamics of classically chaotic Hamiltonians originates in the break- 
down of the classical limit itself for large times. Stated otherwise, the two limiting procedures 
/i — »• and t — *^ cxD do not commute for classically chaotic Hamiltonians P,^JT6[| . How is this 
related to the issue of computability raised above? Long-time unpredictability is precisely 
the condition that characterizes VS as uncomputable, in contrast to VS which is computable 
by way of the Schroedinger equation. Of course, for sufficiently small times the two can be 
arbitrarily close. That a quantum system can temporarily exhibit dynamic instability was 
recognized early in a numerical study of the quantum kicked rotor model [|l7l, and has since 
been investigated extensively [^. 

Indeed it is useful to restate the general arguments above in more explicit terms for this 
widely studied model. Classically, the kicked rotor is a driven, one- dimensional system with 
chaotic orbits for appropriate values of the driving force. The MB orbits for this system, on 
the other hand, are governed by (|T2p, which amounts to a single, non-autonomous first-order 
differential equation. As such, this system cannot possibly be chaotic |1^, implying that 



in the course of time the quantum orbits will deviate arbitrarily far from the classical ones. 
Therefore the classical limit must of necessity fail for this system. Correspondingly, k as 
measured according to Eq. (|13D must vanish. 

In a numerical simulation of the quantum kicked rotor for the purpose of verifying this 
prediction, we integrated Eq. ( p!2|) together with the Schroedinger equation, and measured 
k along the resulting MB orbit using Eq. ([TSp. The result was zero |]rU|, thus verifying 



the predicted regularity of the orbit. An alternative calculation using a pair of nearby 
orbits reconfirmed the regularity of these orbits. This raises the question of whether k can 
ever be positive for this model. To answer this question, we explored the behavior of the 



hybrid quantity /c'^^* = < V^S[q(t, qo, Po),^] >p^,- Observe that here we are considering 



p.v. 



the quantum objects ip and S as evaluated along the classical orbit q(t, qo, Po) (in place 
of MB orbits). Note that while fc'*^'' is not required to vanish by the general arguments 
above, there is no a priori reason to the contrary either, despite the chaotic nature of the 
classical orbit. Numerical results, however, showed nonzero values of k^^'' for a variety of 
conditions, in striking contrast to the behavior of k. This comparison clearly identifies the 
difference between the chaotic classical orbits and their quantum counterparts as the origin 
of the dynamical stability of the quantum kicked rotor. It also indicates that the mean 
decay rate of V^[q(t, qo, Po), ^], the restriction of the wavefunction to the classical orbit, is a 
discriminating footprint of classical chaos in quantum dynamics. A systematic study of this 
intriguing phenomenon must await a better analytical understanding of the hybrid object 

^[q(^,qo,Po),i]- 

Thus far we have considered k for single MB orbits. What is the MB ensemble average 
of k, denoted by k? Recalling that the MB ensemble average is the same as the quantum 
mechanical expectation value, we find from Eq. ( |T5|) or (13) 



limt^^^ J dq\ij{q,t)f\nmq,t)\-\ (16) 



The integral occurring in (^) is recognized as the information associated with the prob- 
abihty density |'?/'(q,t)| . It is actually the entropy associated with an idealized position 
measurement, as shown in Ref. [^. Thus we have arrived at the result that the ensemble 
average of the quantum KS invariant equals the mean growth rate of the position measure- 
ment entropy. Precisely this result was proposed in Ref. |^ on the basis of more intuitive 
considerations. The fact that the original information-theoretical significance of the KS en- 
tropy has naturally re-emerged in the quantum domain may be viewed as an affirmation of 
the present approach. 

As is well known |T7|,^ , short-term measurements of k typically show positive values. In 
light of Eqs. (|T5|) and (p!6|), this suggests a tendency of the probability density to decrease 
along certain MB orbits thereby developing information, i.e., structure and contrast, in 
configuration space. Presumably, points of low concentration would represent the interiors 
of regions of instability corresponding to classical chaos, with the boundaries as points of 
high concentration. This scenario is reminiscent of the scarring phenomena discovered by 



Heller ||2^ and merits further exploration. 
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